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The observed value of the Higgs mass indicates an instability of the Higgs scalar at large energy 
scales, and hence also at large field values. In the context of early nniverse cosmology, this is often 
considered to lead to problems. Here we point out that we can use the instability of the Higgs 
field to generate an Ekpyrotic phase of contraction. In the context of string theory it is possible 
that at very high energy densities extra states become massless, leading to an S-brane which leads 
to the transition between a contracting phase in the past and the current expanding phase. Thus, 
the Higgs field can be used to generate a non-singular bouncing cosmology in which the anisotropy 
problem of usual bouncing scenarios is mitigated. 


I. INTRODUCTION 

Based on the observed value of the Higgs mass [T] , the 
Higgs potential is unstable at large energy scales where 
the quartic self coupling constant and hence the poten¬ 
tial become negative [2], with important consequences 
for cosmology (see e.g. mm- Since large field values 
lead to large energy densities, the instability will arise 
at large field values which are relevant to early universe 
cosmology. It is usually assumed that quantum gravity 
effects will cause the potential to turn positive again at 
very high values (conservatively speaking, at value where 
the Higgs energy density approaches the Planck density). 
There will be a global minimum of the potential which 
corresponds to Anti-de-Sitter space 

This instability is not a problem for the Standard 
Model of particle physics at the present time, since the 
usual Higgs vacua remain local vacua, and their life time 
is much larger than the age of the Universe. However, 
when applied to early universe cosmology some problems 
may arise. In the radiation phase of Standard Big Bang 
cosmology, finite temperature corrections to the poten¬ 
tial are able to stabilize the Higgs field and localize it in 
the usual vacua. The Standard Big Bang Model phase of 
cosmology, on the other hand, must be preceded by an 
early phase which solves the horizon and flatness prob¬ 
lems and allows a causal generation mechanism of fluc¬ 
tuations. The current paradigm for this early phase is 
inflation, an early phase of accelerated expansion [^. In 
simple models of inflation, the energy scale at which in¬ 
flation takes place is of the order of 10^®GeV. Unless 


* Electronic address: rhb@physics.mcgill.ca 
^Electronic address: yifucai@physics.mcgill.ca 
^Electronic address: wanyp@ihep.ac.cn 
^Electronic address: xmzhang@ihep.ac.cn 

^ An easy way to avoid this instability problem is to assume that 
there is new physics at scales lower than the inferred Higgs in¬ 
stability scale which uplifts the potential. In this paper we shall 
assume that this is not the case. 


the Higgs field is directly coupled to the inflaton, the 
field which generates inflation, the Higgs will perform a 
random walk during the period of inflation with a typi¬ 
cal amplitude of H, and time scale of where H is 
the Hubble expansion rate during inflation. This random 
walk will lead to the Higgs field having a finite probabil¬ 
ity of landing in its true Anti-de-Sitter (AdS) vacuum. 
This is a serious problem for our universe [6] (see also 
[7] for earlier discussions of the implications of the Higgs 
instability for inflation, and [S] for related work). 

In this paper we point out that, in the context to alter¬ 
natives to inflation ^ the instability of the Higgs potential 
may be a virtue rather than a problem Specifically, 
we consider a contracting universe described by General 
Relativity coupled to Standard Model matter. Early in 
the contracting phase, the equation of state will be that 
of cold matter, followed by a radiation-dominated phase. 
Eventually, the instability of the Higgs potential will set 
in. As we show, this generates a phase of Ekpyrotic con¬ 
traction, a period which smooths out any pre-existing 
anisotropies. At Planck densities we postulate that a 
new set of fields becomes effectively massless, an effect 
which leads to an “s-brane” in the low energy effective 
action. This s-brane leads to a non-singular transition 
from contraction to expansion, as studied in detail in a 
different context in m- The Higgs field itself bounces off 
a potential barrier and returns back to small field values 
as the universe expands. 

Thus, we show that the instability of the Higgs field al¬ 
lows us to construct a “natural” non-singular cosmology 
which is free from the usual BKL instability [T^] which 
faces bouncing cosmological models m- Since the uni¬ 
verse starts in a matter phase of contraction, our model 
yields a simple realization of the “matter bounce” alter- 


^ See [9] for a review of various alternatives to cosmological infla¬ 
tion. 

^ See also m where the instability of the Higgs was used to con¬ 
struct the contracting phase preceding a bounce. 
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native to inflation [Til ITS] . 


II. MODEL 

We will be considering the Higgs sector of the Standard 
Model of particle physics. The Lagrangian is 

Co = X- V{h ), (1) 

where h denotes the Higgs field, X is the standard kinetic 
term {X = and V{h) is the Higgs potential 

energy. The bare potential energy is given by 

V{h) = \\o{h^-v^)\ (2) 

where Aq is the tree level renormalized coupling constant 
and V is the vacuum expectation value of the field. 

The bare potential is subject to quantum corrections. 
At one loop level, the corrections have been studied in 
detail in m (see also [H])- There are logarithmic cor¬ 
rections to the potentail which enter with a positive sign 
for bosons and with a negative sign for fermions. For the 
observed value of the top quark the fermion contribution 
dominates and the one loop effective potentail can be 
written as 


V^^Hh) = ^X{h){h^-vY, ( 3 ) 


Following m , we assume that at high energy densities 
a new sector of states becomes effectively massless. This 
follows if we consider a superstring model which reaches 
an enhanced symmetry point at some critical density m 
[19]. In this case these states must be included in the 
low energy effective action as a term arising only at a 
particular density, or equivalently at a particular time tc 
when this critical density is achieved. The action S hence 
contains a term of the form 

6S = fi6(t — tc) ■ (6) 

Such a term is called an S-brane, and /i is its tension 
which is set by the mass scale of the new physics. 

An S-brane can be viewed as a relativistic topological 
defect which is space-like. It is well known that for such 
defects the pressure in direction of the defect is negative 
(i.e. positive tension), and the pressure perpendicular to 
the defect vanishes. For an S-brane, this implies that 
the induced energy density vanished and the pressure is 
negative. Hence, for the model we are considering the 
energy density p and pressure p are given by 

p = X + V (7) 

p = X - V - p6{t - tc) ■ (8) 


with 


X{h) 


^0 - Wn(^) , 


( 4 ) 


where 6 is a positive number set by the top quark 
Yukawa coupling constant, and A gives the scale at which 
the instability sets in, which according to the current 
best Higgs and top quark mass measurements is about 
10i°GeV 0. In the following we shall work with the one- 
loop effective potential but we will omit the superscript. 
In our numerical study we will work with a larger value 
of A in order to reduce the magnitude of the hierarchy in 
energy and time scales. 

At field values which correspond to Planck-scale or 
string scale energy densities we assume that there will 
be corrections which uplift the potential. We consider 
(following [4]) the leading order non-renormalizable term 


SV 




( 5 ) 


where the coupling constant g can be absorbed into the 
mass scale M and we can set <7 = 1. In Figure 1 we 
present a sketch of the one loop effective potential in¬ 
cluding the above term. 


From the above it follows that an S-brane leads to the 
violation of the weak energy condition which allows for a 
transition between a contracting phase and an expanding 
phase. 


The cosmological scenario which we have in mind is 
now the following. We begin in a contracting universe 
with the Higgs field close to today’s minimum h = v. 
Initially, h is oscillating about v with an amplitude which 
is small compared to v. The corresponding equation of 
state is that of (when averaged over time) cold matter, 
i.e. p = 0. The amplitude of field oscillations will grow as 
the universe contracts, and eventually it reaches of order 
V. After that point, h is free also to explore negative 
field values. In fact, after some time of contraction the 
local Higgs potential barrier at /i = 0 becomes negligible 
and h will oscillate in a potential which looks quartic. 
At this point the equation of state will change to that of 
(again averaged over time) radiation, i.e. p = The 
early dynamics of h in our model is shown in Fig. 2, as 
follows from numerically solving the equations of motion 
discussed in the following section. 
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FIG. 1: Sketch of the one loop effective potential of the Higgs field including the extra term § coming from postulated 
quantum gravity effects. The instability of the Higgs sets in at field values above A. Note that the local minimum of V(h) 
at h = 0 is not visible on this plot due to the large hierarchy of scales between the electroweak symmetry breaking scale and 
the Planck scale. We have used the following values of the constants appearing in the potential: v = 246GeV, Ao = 0.129, 
b = 0.0187, A = 2.4 X 10^®GeV, <7 = 1 and M = 2.4 x 10^®GeV. The first three of these parameters reflect the measured masses 
of the Higgs and the top quark, the final term reflects the assumption of a quantum gravity-induced wall in the potential at 
values of h corresponding to the reduced Planck mass. 


The amplitude of the (anharmonic) oscillations of h 
will continue to grow until it reaches the local maximum 
of the potential at h ^ A. At that point the instability 
of h sets in and h will start rolling down the potential to¬ 
wards negative values of V. Once the potential becomes 
negative, the equation of state of the Higgs field becomes 
of Ekpyrotic type [20], i.e. 

w = - > 1. (9) 

P 

As was realized in |21] in the context of the matter bounce 
scenario, this phase of Ekpyrotic contraction has the 
virtue of diluting anisotropies [22j . Thus, our scenario 
in a completely natural way solves the main problem of 
a bouncing scenario, namely the anisotropy problem. 

Once the energy density in the contracting phase ap¬ 
proaches the string scale (or Planck scale) density, two 
effects take place. Firstly, the Higgs held hits the “po¬ 
tential wall” where V{h) sharply rises due to the extra 
contribution SV to the potential from ([^. This causes h 
to slow down and V(h) to turn positive. Secondly, at the 
string density the S-brane is encountered. We will 


take this to happen at the time when h comes to rest, 
i.e. X = 0 (we comment later on this assumption). At 
this point, a transition between contraction and expan¬ 
sion takes place. Because for a symmetric bounce the 
kinetic energy is negligibleimmediately before and after 
the bounce the equation of state parameter w will ap¬ 
proach w = —1 from both sides. At the bounce point 
itself, the S-brane leads to a value of w which is formally 
w = —oo (this shows the consistency of our analysis with 
the general theorems of [23] that the single fluid matter 
with a crossing of w = — 1 is required to obtain a bounce 
[ 22 ]). 

After the bounce, the Higgs held h will roll back down 
the potential (towards the origin) picking up kinetic en¬ 
ergy, and then using this kinetic energy to roll back to 
h ~ A and back into the local minimum at h = v. 
The dynamics of h in the phase of Higgs instability and 
around the S-brane bounce is shown in Figs. 3 and 4. 
Once again, these plots are obtained by numerically solv¬ 
ing the equations of motion discussed in the following 
section. 


In the above we have focused on the Higgs sector of the 


Standard Model matter. We should keep in mind that the 
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FIG. 2: The evolution of the Higgs field early in the contracting phase. The initial conditions were taken to he h = v + 5h with 
a small offset 5h = O.ln and h = 0. Initially, h oscillates about h = v, but because of the Hubble antidamping the amplitude of 
oscillations grows and eventually h will be able to cross over the local potential maximum at h = 0. Subsequently, h executes 
oscillations about h = h. Before the transition the time-averaged equation of state is w = 0, afterwards w = 1/3. The figure 
shows the evolution of h (vertical axis) as a function of time (horizontal axis). The field values are in units of v, the time values 
in units of v~^. In this numerical simulation, the h dependence of the quartic coupling constant and the extra h® terms are 
neglected since they are not important during the initial dynamics. Also, in order to see both the increase in the amplitude 
and the oscillation period easily, we have set v = 10“® in Planck units. 


other particles of the Standard Model also contribute to 
the dynamics. If we have in mind starting the contract¬ 
ing phase in a state which looks like the time reverse 
of our present universe, the rest of the Standard Model 
matter will initially be in a cold matter-dominated state 
as well. The radiative degrees of freedom of the Stan¬ 
dard Model matter will lead to a transition between the 
matter-dominated phase and the radiation phase earlier 
than if we only consider the Higgs field dynamics. This 
does not change our scenario. It simply causes the Higgs 
field to climb its potential towards the instability point 

^ A at a different rate. Once the Higgs field starts its 
descent towards negative values, the Ekpyrotic equation 
of state of the Higgs field will cause the Higgs to rapdily 
come to dominate matter, and from then on the analysis 
is exactly how it was described above. 

In the following we will take a closer look at the dy¬ 
namics of the Higgs system. Since in the Ekpyrotic phase 
spatial gradient terms also get washed out we will focus 
on the homogeneous dynamics. 


III. EQUATIONS OF MOTION 


As usual, it is convenient to write the dynamics in 
terms of a rescaled field 


and in terms of conformal time rj which is related to the 
physical time t via 


dt = a{t)dr]. 


( 11 ) 


The Higgs field equation of motion then becomes 

u" - + H')u = -a^Vh , (12) 

where % is the Hubble expansion rate in conformal time, 
namely % = a'/a, and a prime denotes the derivative 
with respect to ry. 

During the radiation phase of contraction the second 
term on the left hand side of (12) vanishes and, in the 


limit \h\^ V the equation reduces to 


u" = -AuM 


(13) 


For field values \h\ <C A, i.e. far from the Higgs instability 
point, we can take A to be a positive constant, and hence 
the solutions yield anharmonic oscillations. 

Beyond the instability point the sign of A changes. As 
mentioned before, the runaway of h to values with neg¬ 
ative potential leads to an Ekpyrotic equation of state 


^ Note that in the contracting phase rj is negative and approaches 

Tj = 0. 


U = ah , 


( 10 ) 
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FIG. 3: The evolution of the Higgs field in the contracting phase after the onset of the instability. The initial conditions were 
taken to be h jnst to the left of the local maximum of the potential with a small initial velocity sufficient to push h over the 
barrier (specifically h = 2 x 10“^ and h — 2 x 10“® in reduced Planck units). The left panel shows the evolution of h (vertical 
axis) as a function of time (horizontal axis). The field and time values in Planck units. To guide the eye we have drawn 
horizontal dashed lines at the value of h corresponding to the onset of the Higgs instability. The numerical result shows that at 
late times the Higgs field becomes localized in the metastable Higgs region near h = 0. The right panel shows the equation of 
state parameter w (solid curve) and the Hubble parameter (dashed curve) (each on the vertical axis) as a function of time. The 
field is seen to rapidly accelerate down the potential, reach negative values of the potential with a resulting Ekpyrotic equation 
of state w ^ 1. After crossing the minimum of the potental, h climbs up the steep potential and comes to rest. At the point 
when h comes to rest we have inserted an S-brane chosen to yield a symmetric bounce, i.e. H simply changes sign, and the 
values of h and h do not change. As can be seen, after the bounce h falls back down to the minimum value of its potential, 
climbs back up towards h — 0, and ends up oscillating about h — 0. At the large field values considered here, the double well 
structure of the potential near h — 0 has a neglible effect. Hence, we have here taken v = 0. 


w ^ 1. This leads to slow contraction 
a(t) - , 


(14) 


with p ^ 1. In this case, the terms of (12) involving "H 


and its derivatives are negligible. During a short time 
interval during which the change in the logarithm in (|^ 
can be neglected the equation of motion takes the form 


// \ 3 

u = Xu 


(15) 


form 

u{r]) = (16) 

V 

with 

f=j- ( 17 ) 

A 

The above rapid growth of h leads to the potential 
energy which scales as 


V 



where A > 0, where has rapidly growing solutions of the 


(18) 
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FIG. 4: This figure is a blowup of the previous one and shows the dynamics of h and w near the bounce point. In particular, 
the figure shows that immediately before and after the actual bounce point w turns negative. The actual numbers show that 
w = —1. The axes and units are as in the previous figure. In this case we do not show the evolution of the Hubble parameter. 


which show that the energy density in the Higgs field 
rapidly comes to dominate over all other forms of en¬ 
ergy density (in particular that of regular radiation which 
scales as a“^). 

The system of equations is completed with the 
Friedmann-Robertson-Walker equations which take on a 
simpler form in terms of physical time. During the stages 
when the Higgs field dominates the energy-momentum 
tensor of matter these equations are 

H^ = ^{X + V) (19) 

where H = a/a is the Hubble expansion rate in terms of 
physical time, and G is Newton’s gravitational constant, 
and 

- = -^{4X-2V- -G)). ( 20 ) 

(X o 


is given by dimensional analysis by 

(23) 

where mg is the string scale and tg is the associated time, 
namely 

t, = . (24) 

Hence 

= G-^'^ml. (25) 

This value has to be compared with the value of H at the 
time when the Higgs energy density reaches the string 
scale density m^. By the first Friedmann equation this is 
given by 

H ~ G'^/'^ml. (26) 


More specifically, inserting the expressions for the equa¬ 
tion for the change in H becomes 

H = -8ttGX + ATTGfi6(t - U ), (21) 


where t(. is the time when the energy density has in¬ 
creased to the point that the S-brane appears. 

From (211 it follows that once the S-brane is hit, the 
value of the Hubble constant jumps by an amount 


AiJ = ^-KG^l. (22) 


In order to obtain a cosmological bounce this jump in 
H has to be large enough to change the sign of H. Let 
us estimate the numbers. For concreteness let us assume 
that the S-brane is due to string states becoming mass¬ 
less, as in the scenario of m- In this case the value of /i 


Thus, we see that the expected jump in H due to the 
S-brane is large enough to change contraction into ex¬ 
pansion. 

Whether the cosmological bounce is symmetric or not 
appears to depend on details of the construction. If the 
time tc is the time when the Higgs field comes to rest 
when running up the potential wall as stringy values, 
then the bounce will be symmetric. If the time tc arises 
earlier, then after the bounce the Higgs field will continue 
to roll up the hill for a while before turning around. Since 
the universe keeps contracting until the S-brane is hit, 
even if h has already turned around, it is also possible 
for the bounce to occur after h has already reached its 
maximal value. 

Let us take a closer look at the matching conditions 
across the S-brane. Given the change of H across the 
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brane by the amount (22), the values of h and h across 
the brane are also determined. We will assume that h 
does not jump. In this case, the equations 


Hi = (27) 

Hi = [F+ + X+] = + XHf , 


where the subscripts — and + indicate the quantities 
right before and right after the bounce. These equations 
determine the value of h after the bounce. For a sym¬ 
metric bounce with X_ = 0 and 77+ = —H_ we obtain 
X+ = 0, i.e. h starts at rest after the bounce. 

In the numerical simulation whose results are presented 
in Figs. 3 and 4 we have solved the equations of motion 
without the brane for times before the bounce is reached. 
We assume a symmetric bounce, i.e. that the bounce 
occurs when h has climbed up the potential at large field 
values to positive values and comes to rest. At this point, 
we reverse the sign of the value of the Hubble parameter, 
as discussed above. After the bounce we again solve the 
equations of motion without the brane. Note that at the 
bounce point w = —1 since A = 0. 

The figures show that the scenario argued for in the 
above approximate analytical considerations is indeed 
obtained. In particular, there is an Ekpyrotic phase of 
contraction followed by a non-singular bounce, and at 
late times in the expanding phase the Higgs field ends up 
in its “regular” vacuum state with \h\ = v. A blowup of 
Figure 4 shows that for the value of A chosen, A = 10“^ 
in Planck units, the Ekpyrotic phase of contraction lasts 
for about 300 Planck times. We have verified numerically 
that the length of the period of Ekpyrotic contraction 
scales roughly as A“^. Thus, for the value A ^ 10“^ in 
Planck units indicated by the current measurements of 
the Higgs and top quark mass, we obtain about 3 x 10® 
Planck times of Ekpyrotic contraction. Whether this is 
sufficient to completely solve the anisotropy problem de¬ 
pends on the initial anisotropy at the beginning of the 
evolution. 

On the other hand, the tuning on the parameter /r re¬ 
quired to obtain a sufficiently symmetric bounce to allow 
the Higgs field to relax to |/i| = u at late times becomes 
more severe the smaller A is. 


IV. COSMOLOGICAL SCENARIO 

The cosmological scenario we have developed is the 
following. We start in a matter-dominated phase of con¬ 
traction with the Higgs field close to its current vacuum 
value V. After some amount of contraction there will be 
a smooth transition to a radiation phase of contaction. 
During both periods h will be oscillating (initially about 
V and later about 0). Once the amplitude of oscillation 
reached a value of order A, the Higgs instability will set 
in. The Higgs field rolls off to negative values of the 
potential, leading to an Ekpyrotic phase of contraction 


during which the energy the Higgs energy density comes 
to dwarf the energy density in all other forms of mat¬ 
ter, and also smooths out anisotropies. Eventually the 
energy density increases to the point when the dynamics 
hits an S-brane, at which point a non-singular transition 
from contraction to expansion sets in. In the expanding 
phase the Higgs field re-traces its evolution in the con¬ 
tracting phase, eventually landing it back in its vacuum 
state h = V. 

This scenario provides a realization of the matter 
bounce alternative to cosmological inflation as a theory 
for the origin of structure in the universe, and as a solu¬ 
tion of the horizon problem. According to this scenario, 
we start early in the contracting phase with fluctuations 
in their quantum vacuum state. The growth of the curva¬ 
ture fluctuations on super-Hubble scales in the contract¬ 
ing phase then transforms the spectrum from a vacuum 
spectrum to a scale-invariant one on length scales which 
exit the Hubble radius during the matter phase of con¬ 
traction [nils]. 

The key question for late time cosmology is whether 
the scale-invariance of the spectrum of curvature fluc¬ 
tuations survives the bounce. This question has been 
studied in many toy models on non-singular bouncing 
cosmology, e.g. in the quintom bounce [Ml ES], in the 
Horava-Lifshitz gravity bonce [55], the bounce [57] ob¬ 
tained in Horava-Lifshitz gravity [55] and in the ghost 
condensate and Galileon bounces [55] (see [501 [QT] for 
the background models), and it was found that on length 
scales larger than the time scale of the bounce phase the 
spectral shape is unchanged. The issue is more subtle, 
however, in the case of an S-brane bounce, as discussed 
in detail in [32] . In this case, the final spectrum depends 
sensitively on the coordinate system in which the brane 
is defined. In the case of purely adiabatic fiuctuations 
the scale-invariance of the spectrum is preserved. 


V. DISCUSSION 

We have shown that the instability of the Higgs po¬ 
tential can have positive consequences for early universe 
cosmology when considered not in the context of the in¬ 
flationary scenario, but in the context of a bouncing sce¬ 
nario. More specifically, we have shown that the insta¬ 
bility of the Higgs field generates an Ekpyrotic phase of 
contraction which smooths out anisotropies and hence 
solves the key problem facing bouncing cosmologies. 

We have argued that an S-brane arising from stringy 
effects can lead to a non-singular bounce. When the S- 
brane is considered to have infinitesimal thickness in tem¬ 
poral direction there will be a jump in the value of the 
Hubble constant, but if the S-brane is smeared out in the 
same way that topological defects are (they have finite 
thickness) then the cosmological evolution is completely 
non-singular. 

Our scenario hence gives a realization of the matter 
bounce alternative to the inflationary scenario of struc- 





ture formation. 

Note that to obtain a symmetric bounce a certain 
amount of tuning of fj, is required. For asymmetric 
bounces the danger is that h will not be able to relax 
to h = ±v at late times in the expanding phase, but 
comes to rest in a AdS minimum. 
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